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Effect  of  the  surface  relaxation  on  spinodal  decomposition  of 
semiconductor  epitaxial  films 

I.  P.  Ipatova , S.  G.  Konnikov,  M.  B.  Lifchits  and  A.  Yu.  Maslov 
Ioffe  Physico-Technical  Institute,  St  Petersburg,  Russia 

Abstract.  Semiconductor  alloys  are  unstable  at  critical  temperatures  because  of  the  positive  en- 
thalpy of  alloy  formation.  The  critical  temperature  of  decomposition  is  the  function  of  elastic 
parameters  which  differ  due  to  structure  relaxation  of  epitaxial  film.  It  opens  the  new  way  to  vary 
the  critical  temperature  and  the  phase  diagram. 


There  are  numerous  experimental  data  [ ] which  show  that  in  the  process  of  semiconduc- 
tor alloy  epitaxial  growth  there  exist  some  temperature  and  composition  intervals  where 
semiconductor  alloy  is  unstable,  and  decompose  into  periodic  structures  with  composition 
modulations.  The  formation  of  such  structures  is  possible  due  to  the  high  value  of  surface 
diffusion  coefficient  (with  respect  to  semiconductor  bulk  diffusion  coefficient). 

To  explain  the  similar  effect  in  the  metal  alloys  Cahn  [ ] and  Khachaturyan  [ ] sug- 
gested the  theory  of  spinodal  decomposition.  Ipatova  et  al  [ ] have  studied  the  spinodal 
decomposition  of  the  thin  epitaxial  films  of  semiconductor  alloys.  It  is  showed  that  for- 
mation of  periodic  structures  is  defined  by  both  the  chemical  energy  and  the  elastic  energy 
of  semiconductor  alloy  caused  by  the  composition  fluctuations.  However  the  phenomena 
of  surface  relaxation  and  reconstruction  which  change  the  crystal  lattice  structure  near  the 
growing  surface  have  not  been  taken  into  account. 

This  paper  deals  with  the  effect  of  the  surface  relaxation  on  the  critical  temperature  of 
spinodal  decomposition.  The  surface  relaxation  means  [ ] that  the  surface  crystalline  plane 
atoms  shift  to  the  new  equilibrium  positions  in  a way  that  crystalline  plane  shifts  as  a whole 
in  the  direction  normal  to  the  surface.  It  leads  to  specific  change  of  the  lattice  parameter 
in  the  direction  of  the  normal.  The  surface  crystalline  planes  of  different  orientation  could 
relax  either  inward  or  outward.  Some  surface  crystalline  planes  do  not  relax  at  all. 

The  semiconductor  alloys  are  unstable  at  some  temperatures  because  they  have  the 
positive  enthalpy  of  formation.  It  means  that  at  T = 0 K the  composition  modulated  alloy 
corresponds  to  smaller  Helmholtz  free  energy  than  homogeneous  alloy. 

Semiconductor  alloys  are  disordered  systems  with  composition  fluctuations  at  each 
point 


c(r)  = c + <5c(r).  (1) 

According  to  the  Vegard’s  law,  the  lattice  parameter  fluctuated  with  composition  (1). 
Coherent  conjugation  of  regions  with  different  lattice  parameters  requires  elastic  strains 
and  results  in  specific  elastic  energy 

h 

^elastic  — ~ J dxdy  J dzkijim  (eij  (r)  - sfj  (r))  (eim  (r)  - s®m  (r))  (2) 

—00 

where  efj  (r)  is  the  static  intrinsic  strains,  A.,y/m  is  the  tensor  of  elastic  moduli. 
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The  elastic  stress  is  equal  by  definition  to 

°ij  (*■)  = J^r}  - (£I>»  (*■)  - Elm  to)  • (3) 

It  follows  from  Eq.  (3)  that  the  stress-free  state  corresponds  to  the  strained  state  of  the 
medium,  sim  (r)  = sfm  (r)  / 0.  The  local  intrinsic  strain  is  a linear  function  of  composition 
fluctuations 

£||  0 0 

e?j(r)  = e?j8c(r),  sfj  — 0 e,  0 (4) 

0 0 £x 

where  ef,  reflects  the  influence  of  relaxation  on  the  symmetry  of  intrinsic  strains.  Disre- 
garding  surface  relaxation:  efj  = eo&ij  ■ 

Semiconductor  alloys  consist  of  binary  compounds  of  similar  chemical  nature.  The 
mixt  ure  components  have  close  chemical  energies.  As  a result,  the  change  of  the  chemical 
energy  in  semiconductor  alloy  is  the  same  order  of  magnitude  as  the  elastic  energy  (2).  The 
total  Helmholtz  free  energy  equals 

8F  = Tjnhom  — Fhom  = SF elastic  T 8 F chem • (5) 


The  requirement  of  total  elastic  energy  minimum  yields  equilibrium  equations  of  the 
theory  of  elasticity,  such  as  the  absence  of  external  forces  in  the  bulk  and  the  boundary 
conditions  of  the  absence  of  external  stresses  at  the  free  surface.  It  is  convenient  to  express 
the  solution  of  these  equilibrium  equations  in  terms  of  the  Green’s  tensor  of  the  semi- 
infinite cubic  elastic  medium  Gap( r,  r')  = Gap{. x —x',  y — y' , z,  z').  This  Green’s  tensor 
was  found  by  Portz  and  Maradudin  [ ].  Finally  we  obtain  the  elastic  free  energy  in  the 
form: 

^Teiastic  — ~ f d3rd3r'Sc(r)  SaaXaappe°^8(r-r')-Aap-^--^—Gap(r,r')  8c(r')  (6) 
Z J L clr a dr p 

where,  considering  the  surface  relaxation: 

Aan  = {cu-cn)1elaelp+cni.c\\-cn){j^s]^(sla+s^)+c\1{^sls)-)8aa8^, 

j <j 

(7) 

(c ii,  c 12,  c'44  being  the  elastic  moduli  in  the  Voigt  notation).  Disregarding  surface  relax- 
ation: Aap  = (Cll  - C12)2£o<W 

Since  elastic  properties  of  semi-infinite  medium  are  homogeneous  in  the  x-  and  _v- 
directions  in  the  substrate  plain,  it  is  convenient  to  use  Fourier  transformation  in  x and  y 
and  express  the  composition  fluctuations  as  a linear  combinations  of  static  concentration 
waves  with  different  kg  = (kx,  ky), 

/i 2 K 

^-^-5c(k|| , z)  exp(/  (kxx  + kyy )) . (8) 

Substitution  (8)  into  (6)  results  in 


NT.14p 


45 


where  S is  the  substrate  area,  the  operator  /ieiastic(kj|,  z,  z!)  is  equal  to 

3 

fieiastic(k||,  2,  z)  = Bq8(z  - z)  + ^ Cs(&)k\\  exp[-avk||  \z  - z |]  + 

5 = 1 

3 3 

+EE^  (®)*H  exp[— q'jA'ii  (h.  - z)]  exp[— ov£||  (h  - z')].  (10) 

,v'=l  ,s=l 

The  first  and  the  second  terms  on  the  right-hand-side  of  Eq.  (10)  correspond  to  the 
operator  B^stic  for  an  infinite  cubic  medium.  The  third  term  is  the  surface  contribution  to 
^elastic-  The  angel  © is  the  angel  between  the  direction  of  the  wave  vector  ky  and  [100]- 
direction  on  the  surface.  The  three  dimensionless  parameters  as  (s  = 1,2,  3)  describe 
the  spatial  attenuation  of  three  static  analogues  of  Rayleigh  waves  generated  by  static 
fluctuations  of  alloy  composition.  These  parameters  are  defined  in  [ ].  The  coefficients 
B( ),  Cs,  Dss/  are  obtained  here  considering  the  surface  relaxation,  and  they  differ  from  the 
ones  corresponding  to  surface  relaxation  disregarding  in  [ ]. 

To  derive  the  criterion  of  absolute  instability  it  is  convenient  to  present  the  variation  of 
the  Helmholtz  free  energy  8F  which  is  quadratic  functional  in  <5c(kj|,  z),  in  the  diagonal 
form.  Thereto  we  define  the  orthonormal  eigenfunctions  cpp (ky , z)  and  eigenvalues  Xp  (ky ) 
of  operator  ^elastic  by  the  following  equation: 

h 

J ^elastic(k||,  Z,  Z )^Pp(k||,  Z )dz  = Xp(]£.\\)(pp (k||,  z)  (11) 

0 

Expanding  an  arbitrary  composition  fluctuation  <5r(r)  in  eigenfunctions  <pp (ky,  z),  one 
obtains  the  variation  8F  of  the  total  Helmholtz  free  energy  of  the  system  in  the  diagonal 
form: 


5F  = l/l^£[(0)+1''(k',]|"''lk')|2-  (12) 

P= 0 

ap(k||)  being  coefficients  in  expansion  of  5c(k||,  z)  in  eigenfunctions  <pp( 'kj,  z). 

It  follows  from  Eq.  (2)  that  the  operator  Elastic  yields  always  a positive  contribution 
to  the  variation  of  free  energy.  It  means  that  all  eigenvalues  Xp  of  the  operator  ^elastic  are 
positive. 

Due  to  mixing  entropy  at  sufficiently  high  temperatures  both  bulk  ternary  III-V  alloys 
and  epitaxial  films  of  alloys  are  stable,  i.e.  ^32//3c2j  > 0,  and  the  chemical  free  energy 
contribution  to  8 F is  positive.  However  the  positive  formation  enthalpies  of  those  alloys 
results  in  the  fact  that  at  sufficiently  low  temperatures  ^32  /'/3e2  j < 0,  and  the  chemical 
free  energy  decreases  due  to  composition  fluctuations.  As  a result,  infinitesimal  fluctuations 
of  composition  in  the  epitaxial  film  may  lead  to  the  decrease  of  the  total  Helmholtz  energy, 
i.e.  to  8F  < 0.  This  means  the  spinodal  instability  in  the  epitaxial  film.  The  negative 

contribution  of  chemical  free  energy,  (d2  f/dc2\  < 0 to  8 F in  Eq.  (12)  competes  with 
the  minimum  of  all  positive  eigenvalues  of  operator  ^elastic  denoted  as  Then  the 
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thermodynamic  criterion  of  the  spinodal  instability  is  as  follows. 


92/ 

3 c2 


+ = 0. 


(13) 


Equation  (13)  determines  the  critical  temperature  of  spinodal  instability  Tc,  the  critical 
alloy  composition  and  the  wave  vector  k||  of  the  "soft  mode"  responsible  for  the  decomposi- 
tion. The  profile  of  composition  fluctuations  in  this  mode  is  described  by  the  corresponding 
eigenfunction  y?o(kj|,  z). 

The  minimum  eigenvalue  A.™1"  of  ^elastic  defining  the  critical  temperature  Tc  became 
different  as  we  consider  the  surface  relaxation.  It  is  shown  that  the  outward  relaxation 
conforms  the  elastic  energy  decrease  thus  the  critical  temperature  of  spinodal  decomposition 
Tc  is  getting  higher;  on  the  contrary  the  inward  relaxation  increase  the  elastic  energy  and 
makes  the  critical  temperature  lower. 

Since  the  effect  of  relaxation  on  the  critical  temperature  of  spinodal  decomposition  is 
different  for  different  surfaces,  it  allows  to  operate  the  miscibility  gaps  and  temperature  Tc . 
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